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ABSTRACT 


Numerical simulations show the formation of self-gravitating primordial disks during the assembly of the first structures 
in the Universe, in particular during the formation of Pop. HI and supermassive stars. Their subsequent evolution is 
expected to be crucial to determine the mass scale of the first cosmological objects, which depends on the temperature of 
the gas and the dominant cooling mechanism. Here, we derive a one-zone framework to explore the chemical evolution 
of such disks and show that viscous heating leads to the collisional dissociation of an initially molecular gas. The effect 
is relevant on scales of 10 AU (1000 AU) for a central mass of 10 Mq (10^ M©) at an accretion rate of 10~^ Mq yr“^, 
and provides a substantial heat input to stabilize the disk. If the gas is initially atomic, it remains atomic during the 
further evolution, and the effect of viscous heating is less significant. The additional thermal support is particularly 
relevant for the formation of very massive objects, such as the progenitors of the first supermassive black holes. The 
stabilizing impact of viscous heating thus alleviates the need for a strong radiation background as a means of keeping 
the gas atomic. 


1. Introduction 


The first stars in the Universe are expected to form from a 
pri mordial gas in minihalos of 10® — 10® at z = 20 — 
30 (lAbel et al.l [200^ iBromm fc LoebI [20031 : lYoshida et al.l 
2008ll . Their formation process is often explored from cos¬ 
mological initial conditions, following the gravitational col¬ 
lapse over many orders of magnitude, to scales of astronom¬ 
ical units (AU) or even below. Due to the characteristic 
temperatures of 300 K in the primordial gas, the first 
stars are expected to be considerably more massive than 
present day stars. Adopting a temperature of ^ 10 K in 
present-day molecular clouds, one may expect an increase 
in mass by a factor of 30^ ® ~ 140 based on simple Jeans 
arguments. 

The actual masses and their distribution depend how¬ 
ever on the details of the star formation process, which are 
only partly understood. A central question concerns the 
role of fragmentation, which may lead to the formation of a 
stellar cluster rather than a central star (IClark et al.ll200§l . 
The latter requires the formation of a self-gravitating disk, 
and its subsequent fragmentation via gravitational insta¬ 
bilities. Indeed, simulations starting from cosmological ini¬ 
tial conditions have now confirmed the formation of self- 
gravi tating disks in minihaloes, as well a s their fragmenta- 
tion ( Claik_et_^2.L 2011 ^ 20121 : iLatif et 

l2013c : iBovino et al.l 2014bl : Susa et al. 2014^ . While these 


simulations followed the evolution for up to 1000 years 
depending on resolution, the Kelvin-Helmholtz timescale 
for a protostar to reach the main sequence is of the order 
of - 10® years, implying that no hnal conclusions can be 
drawn on the properties of the stellar system, and that the 
simulations so far have only explored the evolution of the 
disks at very ear ly stages. An excepti on in this respect 
are th e studies bv iHosokawa et ^ (120111) and iHirano et al.l 
(I2OI4II . who have followed the formation of primordial stars 
until reaching the main sequence, however at the price of 
a 2D-approximation. The latter provides a hint towards 
the final mass distribution, in particular if fragmentation 
is not efficient or if the resulting clumps efficiently merge 
with the central object. For a further discussion concerning 
the state-of-the-art of current simulations expl oring Pop. Il l 
star formation, we refer to the recent review bv iGreij (l2015ll . 


In massive primordial halos of 10^ — 10® M0, the 

5 


has been suggested 

to occur (e.g. Bromm & Loeb 

2003 

; KoushiaoDas et al.l 12004 iLodato & Nataraian 

2006 

; Bea-elman & Sh 

osman 120091: iNataraian 2011; 

Volonteri & BellovarvI 

20121). as potential progenitors 


of the obs e rved s upermassive blac k hole s at z ^ 6 — 7 
(iFan et al.l 20011: Fan et ^ l2004 120061 : iMortlock et al.l 
I2OIII : iBanados et al.l 1201411 . In the presence of a strong 
photo-dissociation background, particular ly high tem- 
peratures of the gas can be achieved (e.g. lOmu 1^ 120011: 
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ISchleicher et alJ|201flh . therefore favoring the formation of 
very massive primordial objects. Simulations have shown 
that self-gr avitating disks form in hal o es cooling via atom ic 


Latif et al. 2013alH : 


Becerra et al.ll201^ . 


hydrogen (iRega ii fc Haehneltl 120091: 

iPrieto et al.l 120131 : iRegan et al.l 120*14 

where central masses of 10 ® Mq can be reached w ithin 
10^ years after the initial collapse (jLatif et al.ll2ni8^ . 

The black hole mass function resulting from such a di¬ 
rect collapse has been derived by iFerrara et al.l (|2014ll us¬ 
ing a semi-analytical approach, considering stellar evolution 
models for rapidly accreti ng protostars l|Hosokawa et al.l 
l2013t ISchleicher et al.ll2013ll . Even for more moderate radi- 
ation fluxes, massi ve stars with 10® —10 ^ Mq may still form 
l|Latif et al.ll20140l . and the rotational support in disks was 
found t o increase for an increa sing amount of molecular 
cooling (iLatif fc Volonteiill2015fl . Also in this context, it is 
therefore important to understand the long-term evolution 
of self-gravitating disks both in the atomic and molecular 
cooling regime. 

To complement such numerical investigations, analyt¬ 
ical studies have been pursued to investigate the prop- 
erties of primordial self-gr avitating disks. For instance, 
iLodato fc NataraianI (l200df) have discussed the role of self- 
gravitational instabilities and fragmentation during the for¬ 
mation of a massive central black hole, which may be ac- 
companied with a starbu rst in the self-gravitating disk. 
Ilnavoshi fc: HaimanI (l2014f) have explored the fragmentation 
properties of a self-gravitating disk cooling via atomic hy¬ 
drogen, concluding that fragments forming via the gravita¬ 
tional insta bilities would subsequently merge with the cen¬ 
tral object. iLatif fc Schleich^ (l2015a^ have extended the 
model into the molecular cooling regime, finding that the 
merging of clumps is likely still efficient, as the characteris¬ 
tic migration timescales of the clumps are shorter than the 
timescales of protostellar contraction. In the regime of very 
massive stars or high accretion rates, the resulting viscous 
heating may even lead to a t ransition from the molecu lar 
to the atomic cooling r egime l|Latif &: Schleich^l2015lJl . 

IFerrara et al.l (|2013l) have explored the stability of self- 
gravitating protogalactic disks in the presence of met¬ 
als, finding that these are subject to strong fragmenta¬ 
tion. The latter may provide an obstacle fo r the forma¬ 
tion o f mass ive black holes, as suggested by iMaver et al.l 
l|2010l l2014f) . In the context of self-gravitating proto¬ 
stellar disks, the role of metal and dust cooling has 
been explored by I Tanaka fc Omul^ (l2014fl . suggesting 
that they may considerably enhance fragmentation. Sim¬ 
ilarly, the results of 3D simulations show that fragmen¬ 
tation is triggered by the additional cooling, especially 
for simulations that are evolved beyond the formation 
of the first p eak. Such beha vi or has been observed 
for ins ta nce bvIClark etliL (l2008ll: I Safranek-Shrade r et al.l 


([201j^ Bovino^t^ 


. , , (2014a): iPeters et al.l l|2014ir and 

ISafranek-Shrader et al.l ( 201511 . 

A substantial debate at this point concerns the ques¬ 
tion of whether a purely primordial and atomic gas is 
required for the formation of supermassive black holes, 
as often assumed in statistical pred ictions for the high- 
redshift black hole popu lation (e.g. iDiikstra et al.l 120081 
l20l4lA'garwal et al.ll20lJl . The typical mechanism consid¬ 
ered here is the photo-dissociation of molecular hydrogen by 
Lyman -Werner radiation, as explored e.g. by IShang et al.l 
([20l0^, which is typically parametrized through a critical 
radiation flux Jcrit- Over the last years, it has been shown 


that the critical flux can be consi derably enhanced by 
considering realistic stellar spectra (ISugimura et al.l|20l4 
lAgarwal et al.l[2015ll . as well as by the dynamics in 3D sim¬ 
ulations, where the initial ionization degree is enhanced by 
shocks and therefore requires an even stronger radiation 
background to prevent the formation of molecular hydrogen 
(jLatif et al.ll2014aLl2ni5all . The simulations suggest a value 
of Jcrit of up to 10®. For the production of the observed 
population of 2 ^ 6 black holes, this value is considerably 
too high, even when considering c hemical uncertainties of 
at most a factor of 5 (|Glovedl2015f) . 

However, it is by no means clear whether such a critical 
radiation flux is in fact required for the formation of m assive 
objects. For instance, iBegelman fc ShlosmanI (|2009ll sug¬ 
gested that massive objects could also form from a molecu¬ 
lar gas in t he presence of self-g ravitating instabilities. Sim¬ 
ulations by ILatif et al.l (l2014hll conhrmed that massive ob¬ 
jects of ~ 10® — 10 '’’ Mq still form for radiation fluxes 
below Jcrit, and accretion rates of ^ 10“® Mq yr“^ can 
be maintained even for a moderate radiation background 
(jLatif fc VolonteIill2015^ . In addition, it is conceivable that 
an ato mic gas forms through alter native pathways. For in¬ 
stance, Ilnavoshi fc Omukail (j2012f) considered the dissocia¬ 
tion of the molecular gas by shocks, and showed that the 
gas w ould then remain a tomi c during the further c o llapse , 
while ISethi et ^ (l2010ll and IVan Borm &: Snaansl (j2013D 
considered a similar effect due to the dissipation of mag¬ 
netic energy. 


More recently, we found indications that viscous heat¬ 
ing in the presence of rotation can substantially heat up 
the interior of a self-gravitating disk, l eading to a transi¬ 
tion to wards an atomic cooling regime (jLatif fc Schleich^ 
l2015bh . The latter can be understood as the viscous heat¬ 
ing rate, which strongly depends on the angular velocity 
of the gas, considerably increases towards smaller scales 
and will ultimately exceed the cooling rate of a molecu¬ 
lar gas. While previous investigations were predominantly 
exploring a fr e e-fall collapse using one-zone model s (e.g . 
Omukai 20011: lOmukai et al.l |2004 I Glover fc Abell l2004 


formation (e.s. Resf 

an & Haehneltl 20091: Glark et al. 2011: 

Greif et al.l 2011L 120121: Latif et al.l l2013dal: iBovino et al. 

2014bl: Prieto et al. 

20131: Regan et al.l 12014; Susa et al. 

2014; iBecerra et al. 

l2015f). we explore here the chemical 


and thermal evolution of such disks at their later stages, 
in particular during the presence of a central massive ob¬ 
ject. For this purpose, we will consider different disk models 
and different chemical initial conditions, representing both 
an initially atomic and an initially molecular gas. For the 
modeling of the chemistry, we employ the p ublicly avail¬ 
able c hemistry package KROMeQ developed by iGrassi et al.l 

(l20Tl . 


The outline of this paper is as follows. In section [2l 
we describe the one-zone model for the evolution in a self- 
gravitating disk in a Lagragian frame as well as the model 
for the chemical evolution. In section [3l we present our 
results for the evolution of disks with an initially molecular 
gas, and section 0] contains the results for initially atomic 
disks. A summary and discussion are presented in sectioiijS) 


^ Webpage krome: http://kromepackage.org/ 
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2. Theoretical framework 

In the following, we outline a one-zone model to describe the 
chemical evolution in the mid plane of primordial disks. For 
this purpose, we assume that the disk is axisymmetric and 
stationary in an Eulerian reference frame. We then consider 
a ring at radius R with mass dM and width dR, implying a 
surface density E = dM/{2TrRdR). In the following, we will 
consider how the gas in the ring moves inward by adopting 
a Lagrangian frame of reference and follow its evolution. 
While the disk is stationary in the Eulerian frame, the gas 
ring will evolve when moving inward, implying that the 
density, temperature and chemical abundances will be time- 
dependent in the Lagrangian frame. 

Depending on the mechanism that provides the effec¬ 
tive viscosity of the disk, we will in the following describe 
two models for the evolution of the gas ring. In our first 
model, we will assume that the viscosity is provided by tur¬ 
bulence and/or magnetic fields, leading to a characteristic 
surface density profile of E cx R~^ in the stationary case, 
as well as a second model which assumes that the viscosity 
is provided by gravitational stresses. The latter requires a 
steeper density profile. 


2.1. Generic disk model 

As outlined above, we consider a radial annulus of mass 
dM and surface density E at radius R in an axisymmetric 
primordial disk assumed to be stationary in the Eulerian 
frame. As a result, the accretion rate M is independent of 
time and position within the disk. The sound speed in the 
annulus is given as 


Cs = ^/^ksT/m^ ( 1 ) 


with ks the Boltzmann constant, T the gas temperature 
and m the mean molecular mass. In the following, we need 
to construct a dynamical model for the evolution of the 
annulus, in particular the radius R and the surface density 
E. The latter is necessarily approximate, as we are not 
considering a fully Keplerian disk, but rather a situation 
where the disk is also stabilized by turbulent pressure, as 
reflected in the rather thick disks forming in primordial 


environments (see e.g. Turke^_^ 20091 : 1 Clark et al.l 120111 : 


ICreif et al.ll2012l : iLatif et al.l2018c l. 


In the generic disk model, we will initially consider a 
situation with E oc R~^, as in the case of a so-called Mestel 
disk. The latter corresponds to a steady-state solution for 
a strongly self-gravitating disk with con stant sound sp eed 
Cs and a constant Toomre-Q parameter (lToomrelll964ll 


Q 


CgGl 

ttGE 


( 2 ) 


of the order 1 (iBertinlflOOTl : lLodatoll200i11 . Strictly speak¬ 
ing, a Mestel disk further extends to infinity, while we are 
here considering disks of a finite size. The latter implies 
that also the rotation curve derived from the Mestel disk 
holds only in an approximate way. While the assumptions 
employed here are certainly an approximation, for instance 
the surface density distributio n in realistic disks is often 
similar to the R~^ profile fe.g. iBeckwith et al.lll99dll . and 
we therefore employ it here in our approximate model. The 


evolution equation of the disk surface densit y in an axisym¬ 
metric model is now given as (|Lodatdl2007f) 




1 9 [ 1 d 

Rm 


(i/Ei^^D') , 


(3) 


where the prime ' denotes the derivative with respect to R. 
Under the assumption that v and U are independent of E, 
which certainly is an idealization, it can be shown that the 
disk surface density evolves on a timescale of the order of 
the viscous timescale, which is given as 

tvis — ^ : ( 4 ) 


where zz the effective viscosity of the gas. Considering the 
disk annulus in our model, the evolution of the surface den¬ 
sity is thus approximately given as 



(5) 


In the case of a Keplerian disk, this expression becomes 
accurate by including an additional factor of 3/2. We will 
however consider more generic situations, where the disk 
is also stabilized by turbulence, and not necessarily domi¬ 
nated by a central source. As we assume here a disk profile 
with E oc R~^, the radius of the annulus must follow a 
similar evolution equation given as 

R = -^. (6) 

^vis 


For a stationa ry thin disk, the continuity equation becomes 
(lLodatoll2007^ 


M = 


dlnU 

din A 


2'KiyT,. 


(7) 


For a generic disk model, the factor |(ilnU/dlni?| will usu¬ 
ally be of order 1, and becomes equal to 3/2 for a Keplerian 
thin disk. We therefore employ the approximate expression 

M = Stti/E. (8 ) 


In our model, we assume that the constant accretion rate 
M is known, and the evolution of the disk surface density 
follows from Eq. O We can therefore employ Eq. |S]to solve 
for the disk viscosity required to maintain the accretion 
rate, yielding 


M 


(9) 


While this identity follows from the assumption of sta- 
tionarity and axisymmetry, such conditions can only be 
achieved if a physical mechanism is present to provide the 
effective viscosity v. We note that the formulation em¬ 
ployed here is clearly approximate, and the evolution in 
the disk can depend both on the source of the viscosity, 
whi ch can be due to turbulence and magnet ic fields (see 
e.g. iBalbiis fc PaDaloizodll999l: iHawlevI 20fl(lll or the pres¬ 
ence of gravitational instabilities (|Toomrdll96^ . We as¬ 
sume here that the Toomre Q parameter is at least initially 
of order 1, implying that self-gravity is relevant when the 
annulus forms. If the gravitational force is balanced by the 
centrifugal force, the Keplerian angular velocity in a disk 
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dominated by self-gravity is given as Q,k = \J‘I'kGYj jR 
(|Lodatd l2007f) . In a realistic disk, we however expect it 
to be partly supported by thermal and turbulent pres- 
sure, as also dem onstrated via numerical simulations (e.g. 
iLatif et ^l2013aO ). For the disk-dominated case, we there¬ 
fore adopt 


I^disk — 


27rGS 


R 


( 10 ) 


where describes the deviation from a fully rotation- 
ally supported disk and G denotes the gravitational con¬ 
stant. We will in the following assume a generic value 
tK = 50%, as primordial dis ks typically have s cale heights 
with g/fl betwee n 0.3 — 1% (iTu rk et al1 |2009l : IClark et al.l 
l2011t[&eif et al.ll2012l : iLatif et al.ll2013al lcl. As we assume 
that S oc the latter implies that the Toomre Q pa¬ 
rameter remains ^ 1 if the sound speed Cs remains approx¬ 
imately constant. We note that our surface density profile 
further implies a non-zero enclosed mass in the innermost 
region of our disk. For the disk models considered here, 
the maximum mass within the central 0.1 AU is 0.003 Mq, 
which is expected not to be relevant for the dynamical evo¬ 
lution. 

If the disk is dominated by the central source, we simi¬ 
larly assume that the angular velocity is given as a fraction 
tK of the Keplerian rotation, implying that 




source 



( 11 ) 


with M* the mass of the central source. In this case, it 
follows that Q oc oc R~^/'^/R~^ oc i.e. Q is 

increasing with decreasing radius, and the interior of the 
disk becomes gravitationally stable. As long as a sufficient 
viscosity is provided by turbulence and magnetic fields, the 
latter is not a problem for the assumed stationarity of the 
disk. In the next subsection, we will consider the possibility 
that the viscosity is provided by self-gravitational instabil¬ 
ities, which will require a steeper relation between S and 
R. We now calculate the disk height as 

^ = (12) 

The latter expression hol ds for gravita tionally stable or 
marginally unstable disks (|Lodatoll2007fl . and is thus valid 
both for Q \ and Q > 1. The number density in the mid 
plane of the disk follows as 


S 

“ 2Hm" 


(13) 


and the mass density is given as p = nm. With these 
quantities, we can als o evaluate the viscous heating ra te. 
The latter is given as (|Lodatoll2007l : iFerrara et al.l[20T^ 


Q+ = iyT,{Rn')^. 


(14) 


The viscous heating rate therefore generally depends on 
the rotational profile. In case of rotation around a central 
source, the viscous heating rate can be evaluated as 


Q + ,source — ^^^^source' (1®) 
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For a disk dominated by self-gravity, on the other hand, we 
have n oc i/S/i? oc R~^, implying that 

Q+,disk = *^S^disk- (1®) 


In this framework, we do not impose thermal equilibrium, 
but rather follow the non-equilibrum evolution of the gas 
temperature in our annulus using the detailed chemical 
model described in section 12.31 The corresp onding evo¬ 
lution equation is given as (jGrassi et al.ll201^ 


dt 


(7-1) 


r-A 

fcsn ’ 


(17) 


where n is the number density of the gas following from disk 
surface density, disk height and chemical composition, F is 
the total heating rate including viscous heating, and A the 
total cooling rate. A further description of these expressions 
is given in section [531 While following the non-equilibrium 
evolution of the system is in principle more accurate than 
the assumption of thermal equilibrium, we note that the 
resulting state generally leads to a situation close to thermal 
equilibrium, i.e. with F ~ A. 


2.2. Self-regulated disk model 

As already shown above, a disk dominated by a central 
source is gravitationally stable in the interior if the evolu¬ 
tion of the annulus is dictated by Eqs. 0 and [HI implying 
E oc R~^. Such a scenario can thus only be maintained if a 
sufficient viscosity is provided by turbulence and magnetic 
fields, as assumed in the generic disk model. However, it 
is conceivable that gravitational instabilities are required 
to transport the angular momentum, and they were also 
shown in pr evious studies to pro vide the stronger contri¬ 
bution fe.g. iFromang et al.l 120041) . While the reality may 
lie somewhere in between, we will explore here the extreme 
case where self-gravity is necessary to provide the viscosity, 
requiring that Q = 1. 

In this model, we assume that the disk is gravitationally 
unstable on scales larger than the current radius R, imply¬ 
ing that gas is transported onto the annulus from larger 
scales, and the surface density of the annulus keeps increas¬ 
ing as in Eq. [SI The radius R however remains constant as 
long as Q > 1, implying R = 0, and only evolves as in 
Eq. [5] if <5 < 1. This condition will effectively alter the 
relation between E and R, and increase the surface density 
until a state of marginal stability is reached. As a result, 
one obtains a higher surface density and number density 
on a given scale, and the ability of the gas to cool increases 
compared to the generic disk model outlined above. 

2.3. Chemical model 

To study the chemical evolution for the disk models out¬ 
lined above, we employ the publicly avai l able c hemistry 
package KROME developed bv lOrassi et al.l ([203). In this 
fra mework, we adopt a network similar to the one described 
bv ILatif et al.l (|2015all . including the species H, H+, H”, 
HJ, H 2 , He, He"*", He"*"^ and e“. The photo-rates are dis¬ 
carded here, due to the efficiency of self-shielding at high 
densities. The network includes the currently most accu- 
rate ra t e coeffi cient for the 3-body H 2 formation derived by 
iForrevI l|2013f) from o n a quantum-mechan ical calculation 
(see also discussion bv iBovino et al.l[2r)14cll . 
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Lenzuni et al. 


Grassi et al.l 

2014 . 


The cooling functions in KROME include the atomic line 
cooling, recombinati on c o oling and Bremsstrahlung cool¬ 
ing as des cribed by Cen (Il992ll . the H 2 cooling function 
derived bv iGlover fc Abell ( 2fl08ll as well as the Ho forma - 
tion heating and cooling as described by lOmul^ (|2000t) . 
Continu u m coo ling processes are treated as described by 
0mukai|jl2000() with the continuum opacities derived by 
(ll99lH . em ploying a new fit pr ovided by 
The iLenzuni et alJ (|l991h opacities 
include all relevant continuum processes, specifically the 
bound-free absorption by H and H“, free-free absorption 
by H, H“, H 2 , H^, H^, H;^, He and He“, Rayleigh scatter¬ 
ing by H, H 2 and He and collisionally-induced absorption 
by H 2 . The viscous heating is included as described in sec¬ 
tion!^ 

The chemical model described here will be called within 
each dynamical timestep of the disk model outlined above 
to follow the chemical and thermal evolution. As a result, 
we obtain the evolution of the chemistry while the gas is 
moving inward. The latter corresponds to a time sequence 
in the Lagrangian frame of the annulus, but also describes 
the structure of the stationary disk in the Eulerian frame. 
In the following, we will typically show the results as a 
function of radius to illustrate the resulting structure of 
the disk. 


3. Chemical evolution for an initially molecular gas 

In the following, we present the results for the chemical 
evolution in primordial viscous disks for the scenarios out¬ 
lined above. For the chemical initial conditions, we consider 
an initially molecular gas. We adopt a helium fraction of 
7.75% with respect to the total number of nuclei. The typ¬ 
ical densities at the beginning of the calculation are of the 
order ~ 10® — 10^^ cm“^. For a conservative assessment for 
the impact of viscous heating, we assume here that the gas 
is very close to fully molecular, with an atomic hydrogen 
fraction of only 10“®. For the electron/proton abundance, 
we adopt a gener ic value of 10~®, w hich can be expected 
at such densities (|Bovino et al.l[2M^ . The further species 
abundances like H~ and Hj are initially set to zero, but al¬ 
most instantaneously reach their equilibrium values during 
the chemical evolution. We also checked that the results do 
not strongly depend on these assumptions. 

In the following, we will present the chemical evolution 
for such an initially molecular gas for the case of a generic 
disk model, both within and without a central source, as 
well as for a self-regulated disk model, for which we pre¬ 
dominantly focus on the case with a central source. 

3.1. Results for a generic disk model with no central source 


Model 

M [Mo yr-i] 

E [g cm 

R [AU] 

T\K] 

NSO 

If?’ 

100 

1000 

300 

NSl 

10-1 

50 

1000 

300 

NS2 

10-® 

25 

1000 

300 

NS3 

10-3 

12.5 

1000 

300 


Table 1. Generic disk models with no central source and an 
initially molecular gas. 


As already shown in section 12.11 the surface density 
scales as R~^ under these conditions, and the surface den¬ 
sities reach peak values between 10® and 10® g cm“® at 
R = 1 AU. The results for this case are given in Fig. [TJ 
We find that the temperature generally increases towards 
the interior as a result of viscous heating <5+ oc 
(top panel). As oc (Eq. [5]), the latter essentially 
scales as oc E/i? oc R~^, implying a steep increase 

as a function of radius, which is only partly compensated 
by the increasing number density of the gas. In the regime 
considered here, we note that H 2 is already thermalized, 
so that the cooling rate only scales linearly with uhsi 
optically depth effects are increasing with density. The 
gas temperature therefore increases towards values between 
1000 — 2000 K, and then rises more gradually. 

This behavior is generally found in all cases consid¬ 
ered here, with minor deviations for an accretion rate of 
0.001 Mq yr“^, where the gas temperature slightly de¬ 
creases again near 0.1 AU, as viscous heating is reduced 
in this case and the collisionally-induced emission (CIE) of 
molecular hydrogen simultaneously becomes relevant. We 
further note that the temperature of the gas increases with 
the accretion rate. This is due to the proportional increase 
in viscous heating, which is only partly compensated by 
the increased number densities and the ability of the gas to 
cool. The viscous heating has therefore a relevant impact 
on the thermal evolution of the gas, even though the dif¬ 
ferences in the temperature are here smaller than a factor 
of 2. 

The chemical abundances for a characteristic case with 
an accretion rate of 0.001 Mq yr“^ (model NS3) are shown 
in Fig. [T] (mid panel). We find that the gas is fully molecular 
during the entire evolution. The atomic hydrogen abun¬ 
dance is initially constant, temporarily increases towards 
10“® around R ^ 1 AU as a result of the temperature in¬ 
crease, and decreases again in the interior, reflecting the 
slightly decreasing temperature. Similarly, the proton and 
electron abundance decreases towards the interior as a re¬ 
sult of the increasing densities. 


For a generic disk model with no central source, we explore 
a range of scenarios with different accretion rates, varying 
between 10~® — 10® vr~^, ranging from typical Pop. HI 

star formation (I Abel et al.ll2002| : Bromin fc : Loeh| 20031) to 
the conditions in atomic cooling haloes ( Latif et al.ll2013(l 
iRegan et al.ll20l3l . We assume the initial annulus to be at 
a radius R = 1000 AU, with a characteristic temperature of 
^ 300 K. The su rface densities of the disk a re determined 
via the model of iLatif fc Schleiched (l2015bll . implying an 
initial condition with Q = 1. The details of the models are 
given in Table [TJ 


The resulting contributions to heating and cooling are 
given in Fig. [1] (bottom panel), again for the case of NS3. 
We note that the viscous heating rate is almost balanced by 
the molecular hydrogen cooling throughout most of the evo¬ 
lution between 1000 and 1 AU, and subsequently the con¬ 
tinuum cooling dominates in the interior. Within ~ 1 AU, 
we also find strong chemical heating and cooling, which 
are however tightly balanced, as the chemistry is close to 
its equilibrium value. During most of the evolution, the 
disk remains in the molecular cooling regime, without any 
strong transitions in the gas temperature. 
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3.2. Results for a generic disk model with central source 

Now, we consider the case of a generic disk model with 
a central source. For this purpose, we consider stellar 
masses from 10 — 10^ Mq, ranging from typical Pop. Ill 
to supermassive primordial stars, and accretion rates of 
10“^ — 10° Mq yr“^. We start following the evolution 
at i? = 1000 AU, taking the initial s urface density and 
gas tem perature from the disk model of iLatif fc Schleicheil 
l|2015lJ) . The details of the models are given in Table [21 


Model 

M, [Mq] 

M [Mq yr“°] 

E [g cm“2] 

r[K] 

SlAl 

10 ^ 

10 “^ 

10 

300 

S1A3 

10 ° 

10 “° 

10 

300 

S2A1 

10 ^ 

10 “° 

40 

500 

S2A3 

10 ^ 

10 “° 

40 

500 

S4A0 

10 ^ 

10 ° 

790 

1900 

S4A1 

10 ^ 

10 “° 

630 

1100 


Table 2. Models for a generic disk model with central source 
and an initially molecular gas. 

In this case, the surface density still scales as R~^, 
reaching peak values of 10° — 10^ g cm“^. The temperature 
however changes more significantly during the evolution. 
For a characteristic case with a stellar mass of 10 Mq and 
an accretion rate of 10“° Mq yr“° (model SlAl), the tem¬ 
perature increases from 300 K at i? = 1000 AU to ~ 2000 K 
at 10 AU, and then raises rapidly to ^ 7000 K after only 
a minor change in radius as shown in Fig. |2] (top panel). 
This behavior reflects the steep dependence of the viscous 
heating rate with radius, scaling as oc The 

scaling relation is thus considerably steeper than for the 
case without a central source, and increases further with 
increasing stellar mass. Here, the temperature increases 
towards a value of ^ 2000 K, where the collisional dissocia¬ 
tion of H 2 becomes relevant. We note that the latter has an 
exponential dependence on the temperature, and therefore 
the molecular hydrogen becomes fully dissociated after the 
temperature increases further. At that point, the gas is no 
longer able to cool via molecular hydrogen, and therefore 
increases until the heating is balanced via atomic cooling. 
The temperature then remains high at 6000 — 8000 K on 
scales below 10 AU. 

In general, the evolution is rather similar also for the 
other cases considered. We note that the transition towards 
the atomic cooling regime occurs earlier for higher stellar 
masses and/or higher accretion rates, as both lead to an in¬ 
crease in the amount of viscous heating. For a 10^ Mq star 
with an accretion rate of 1 Mq yr“°, the viscous heating 
is particularly efficient, leading to an almost instantaneous 
transition already at i? ^ 3000 AU. This transition to¬ 
wards the atomic cooling regime may considerably increase 
the stability of the accretion disk in t he environment of 
superm assive stars, as also described bv iLatif fc Schleiched 

||20T^ . 

The chemical abundances for a representative case 
with a stellar mass of 100 Mq and an accretion rate of 
10“° Mq yr“° (model S2A3) is given in Fig. |2] (mid panel), 
clearly showing the transition from the molecular to the 
atomic regime at i? ^ 30 AU. The latter is accompanied by 
a significant increase of the ionization fraction due to the 
higher temperature. We note that the molecular hydrogen 


abundance starts increasing again with density after the 
transition, but does not reach values where the H 2 cooling 
would be significant. In Fig. |2] (bottom panel), the heating 
and cooling rates are plotted for this reference case, con¬ 
firming that the cooling is initially dominated by molecular 
hydrogen, and the atomic line cooling takes over after the 
transition. Within the central 0.3 AU, we note that also 
chemical heating due to H 2 formation starts getting signif¬ 
icant, but is still balanced by the atomic hydrogen cooling. 
The central temperature then raises to a value of ^ 9000 K 
at 0.1 AU. 

3.3. Results for self-regulated disks 

As already discussed in section 12.21 it is conceivable that 
gravitational instabilities are required to provide a suffi¬ 
ciently large effective viscosity to transport the angular mo¬ 
mentum, in particular in the regime of high accretion rates. 
We therefore consider the evolution in a self-regulated disk, 
starting with an initially molecular gas. We predominantly 
focus on the case with a central source. As we mentioned 
in section the evolution in the case without a central 
source would be very similar to the generic disk model. We 
explore characteristic masses of the source of 1 — 10° Mq, 
including typical Pop. Ill stars and supermassive stars. Of 
course, even stronger effects can be expected in case of more 
massive stars. 

We focus here on typical accretion rates of 
10“° Mq yr“°, which is particularly relevant for the 
formation of supermassive stars in atomic cooling halos, 
and where the heating is particularly strong. We follow the 
evolution of the disk from a radius of i? = 1000 AU, and 
evaluate the initial surface density an d gas te mperature 
using the model of iLatif fc Schleich^ (l2015bll . Due to 
the generally higher surface and number densities in this 
scenario, we restrict the evolution here to a minimum 
radius of 1 AU. In fact, considering stellar evolution 
models for supermassive stars, their typical ex tend may 
even exceed such a scale (iHosokawa et al.l 120131) . Overall, 
a summary of the scenarios considered here is given in 
Table 01 


Model 

M, [Mq] 

M [Mq yr“°] 

E [g cm“2] 

T[K[ 

SI 

1 

10 “° 

2.5 

300 

S2 

10 

10 “° 

10 

400 

S3 

10 ^ 

10 “° 

40 

500 

S4 

10 ° 

10 “° 

160 

800 


Table 4. Models for a self-regulated disk with initially molecular 
gas. 


For the self-gravitating disk models, due to the re¬ 
quirement of Q = 1, the resulting relation between sur¬ 
face density S and radius R is considerably steeper than 
for a generic disk model, requiring roughly a relation of 
E oc i?“° °. From the results given in Fig. 0] (top panel), 
it is evident that characteristic features occur when the 
temperature increases strongly, i.e. at the transition from 
the molecular to the atomic cooling regime. In this case, 
the higher gas temperature needs to be compensated with 
a higher gas surface density to ensure the condition that 
Q = 1, thus temporarily steepening the relation between 
surface density and radius. 
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Abbreviation 

cooling and heating terms 

reference 

viscous 

viscous heating 

Lodato 120071 

chem. 

net chemical heating / cooling due 

Omukai 120001 

chem. heat 

chemical heating due to H 2 formation 

Omukai 120001 

chem. cool 

chemical cooling due to H 2 dissociation 

Omukai 120001 

atomic 

atomic line cooling, recombination cooling, Bremsstrahlung 

Cen 119921 

H 2 

H 2 line cooling 

Glover & Abel 120081 

cont 

continuum cooling 

Lenzuni et al. 119911 


Table 3. Cooling and heating functions listed under different abbreviations in Figs. [T] [313111] and [S] 


For a characteristic case with 1 Mq and an accretion 
rate of 10“^ Mq yr“^ (model SI), the temperature in¬ 
creases first from 400 K at 1000 AU to ^ 1800 K at 30 AU, 
and subsequently increases more gradually, due to the col- 
lisional dissociation cooling of molecular hydrogen, towards 
temperatures above 2000 K near 3 AU, where H 2 becomes 
fully dissociated and the temperature increases to ~ 5000 K 
(Fig. [21 2nd panel). A similar behavior is shown for higher 
stellar masses. This is however more extreme, implying that 
the transition to the atomic regime occurs earlier, with a 
higher characteristic temperature. For instance in the case 
of a supermassive star with 10^ Mq, the transition to the 
atomic regime occurs already at 300 AU, with a peak tem¬ 
perature of 7000 K. 

For our reference model SI, the evolution of the species 
abundances is given in Fig. [3] (3rd panel), confirming the 
transition from the molecular to the atomic regime near 
10 AU as well as the subsequent build-up of molecular 
hydrogen due to the high densities of the gas. It there¬ 
fore starts affecting the cooling again around densities of 
10^® cm“®. The viscous heating is initially balanced by the 
molecular cooling between 30 and 1000 AU (Fig. [21 bottom 
panel). Between 30 and 3 AU, the main cooling channel 
is due to H 2 collisional dissociation cooling. In the central 
3 AU, the contributions of viscous heating, continuum cool¬ 
ing as well as chemical heating and cooling strongly balance 
each other, and molecular hydrogen line cooling starts play¬ 
ing some role due to the increasing H 2 abundance. How¬ 
ever, it is also visible for instance in the case S2 that the H 2 
abundance is collisionally dissociated, as the critical tem¬ 
perature for H 2 dissociation decreases with density and the 
disk evolves into a denser regime. We expect similar effects 
to occur in the scenarios S3 and S4, where the calculation 
becomes however numerically unstable at this point as a 
result of the stronger viscous heating. 

4. Chemical evolution for an initially atomic gas 

Now, we consider a case where the gas in the disk is initially 
atomic. In this case, we adopt an initial molecular hydrogen 
abundance of 10“® and an initial temperature of 10^ K. Due 
to the higher temperature, we also assume an increased 
ionization degree of 10“®. In the following, we present the 
chemical evolution for an initially atomic gas both in the 
case of a generic disk model, as well as for a self-regulated 
disk. 

4.1. Results for a generic disk model with and without 
central source 

For the case of a generic disk model, we consider scenar¬ 
ios with and without a central source, with source masses 


of 10 — 10"^ Mq. We adopt here an accretion rate of 
10“^ Mq yr“^ to explore a case where viscous heating is 
particularly relevant. We start with a characteristic surface 
density of 10 g cm“^ at 1000 AU. We note that these disks 
are not necessarily self-gravitating, in particular in the pres¬ 
ence of a central source, and the viscosity to build up the 
disk thus needs to be provided by turbulence or (magneto- 
)hyrodynamical instabilities. A summary of the explored 
models is given in Table [SJ 


Model 

M, [MqI 

M [Mq yr-i] 

S [g cm 2] 

T[K] 

HO 

0 

10 -^ 

10 

101 

HI 

10 

10-1 

10 

lOi 

H2 

100 

10-1 

10 

lOi 

H3 

10 ^ 

10-1 

10 

lOi 


Table 5. Models for a generic disk model with an initially 
atomic gas. 


The evolution of the gas temperature as a function 
of radius is given in Fig. [3] (top panel). The tempera¬ 
ture remains high throughout the entire evolution, from 
R = 1000 AU to R = 0.1 AU, but decreases slightly from 
initially 8000 — 9000 K to 6000 — 9000 K in the interior. We 
note that the tem-perature slightly increases with central 
mass as a result of the enhanced viscous heating. The role of 
viscous heating is overall less pronounced co mpared to the 
i nitial ly molecular regime, as also found by [Ferrara et al.l 
(l2013f) . In fact, the temperature remains high even in the 
absence of a central source, though may somewhat increase 
in the presence of viscous heating. 

For a characteristic case with a central source of 10 Mq 
and an accretion rate of 10“^ Mq yr“^ (model HI), the re¬ 
sulting species abundances are plotted in Fig. [3] (mid panel). 
We note that the gas is atomic in the entire regime, and 
H 2 remains strongly suppressed via collisional dissociation 
due to the high temperature. Within the central 10 AU, 
the H 2 abundance starts increasing due to a mild decrease 
of the temperature, however leading only to an abundance 
of lO”"* at 0.1 AU, which is insufficient to drive the cool¬ 
ing. The ionization degree only evolves very gradually as 
a function of radius, and is slightly decreasing due to the 
increasing densities. The viscous heating is balanced by the 
atomic cooling through almost the entire evolution (Fig. (H 
bottom panel). Only in the very central 0.3 AU, chemical 
heating via three-body H 2 formation becomes important as 
well, but is still balanced by the atomic line cooling. The 
molecular hydrogen cooling remains negligible throughout 
the entire evolution. 


Article number, page 7 of ll3l 





















4.2. Results for self-regulated disks 

Now, we finally explore the chemical evolution of initially 
atomic disks in a self-regulated scenario with Q = 1. In this 
regime, we focus again on the case with a central source, as 
otherwise the evolution would be very similar to the case of 
a generic disk model. We consider here central sources with 
masses of 1 — 10 ^ Mq, and characteristic accretion rates of 
10“^ Mq yr“^. We start with an initial surface density of 
10 g cm“^, which however adjusts early in the evolution 
due to the requirement of Q = 1. Due to the generally 
higher surface and number densities in this scenario, we 
restrict the evolution here to a minimum radius of 1 AU. In 
fact, considering stellar evolution models for supermassive 
stars, their typical ex tend may even exceed such a scale 
(jHosokawa et al.ll201,‘lll . The overall models considered here 
are summarized in Table [B1 


Model 

M, [Mq] 

M [Mq yr-1] 

S [g cm 

HSO 

1 

10-^ 

10 

HSl 

10 

10-1 

10 

HS2 

100 

10-1 

10 

HS3 

10^ 

10-1 

10 


Table 6. Models for a self-regulated disk with initially atomic 
gas. 


The surface density of the gas as a function of radius is 
given in Fig.(top panel), showing that the surface density 
profile clearly steepens compared to a generic disk model, 
requiring a scaling as E oc Different from the case 

with initially molecular cooling, the profile strongly resem¬ 
bles a power-law with no particular features, as the tem¬ 
perature evolution is generally more gradual and mostly 
without strong transitions. The temperature initially re¬ 
mains high, while some continuum cooling may set in be¬ 
tween 60 and 10 AU (Fig.jSJ 2nd panel). The temperature 
then decreases more steeply towards ~ 3000 K. For central 
sources with 1 and 10 Mq, we find that the temperature 
subsequently increases again, as a result of viscous heating 
and II 2 collisional dissociation. For higher central masses, 
the calculation becomes unstable at the point of this transi¬ 
tion, due to the correspondingly higher heating and cooling 
rates. In principle, we expect however a similar evolution. 

The chemical evolution is given for a characteristic case 
with a central source of 1 Mq and an accretion rate of 
10“^ Mg yr“^, corresponding to model HSO (Fig. [5l 3rd 
panel). While the gas is atomic in the entire range, the H 2 
abundance increases significantly from ^ 10“® at ^ 50 AU, 
reaching a peak value close to 0.1 at 1 AU. The ionization 
degree initially decreases only slightly with density between 
1000 and 10 AU, and then drops more significantly towards 
10 “^^ due to the lower temperatures. 

As shown in Fig. [S] (bottom panel), the viscous heat¬ 
ing is balanced by atomic hydrogen cooling between 1000 
and 30 AU. In the interior, the contribution of the contin¬ 
uum cooling is subsequently close to the viscous heating. 
In this regime, the chemical heating and cooling as well as 
the contributions of viscous heating and continuum cooling 
strongly balance each other, and the temperature decrease 
is driven by the molecular hydrogen cooling. The charac¬ 
teristic temperature for H 2 dissociation however depends 
on the density, and the critical point for H 2 dissociation 


is indeed reached at i? ^ 1 AU, leading to the dissocia¬ 
tion of molecular hydrogen and a resulting increase of the 
temperature. 

Also in this case, the role of viscous heating is therefore 
less pronounced for the initially atomic than for the initially 
molecular regime. However, due to the higher gas densities 
in the self-gravitating regime, the H 2 formation is enhanced 
for these cases, so that the cooling is initially somewhat en¬ 
hanced, and the viscous heating prevents the gas from be¬ 
coming molecular again in the central 1 — 10 AU, therefore 
somewhat contributing to the stability in the very central 
region. This behavior concerns relatively small scales. A 
supermassive star with 10^ Min may already have radii 
up to ^ 100 AU (|Hosokawa et al.ll2013^ . providing a nat¬ 
ural cut-off for the inner disk. It thus appears likely that 
the viscous heating has only a minor effect in the initially 
atomic regime. 


5. Summary and discussion 


We have developed a one-zone framework to describe the 
chemical and thermal evolution of primordial disks depend¬ 
ing on the mechanism which provides the viscosity of the 
disk. Assuming a sufficient viscosity is always available 
as a result of turbulence and/or magne tic fields (see e.g. 
iBalbus fc PaDaloizoiJll99^lHawlevll200flh . we have derived 
the generic disk model under the assumption of station- 
arity and axisymmetry as well as a surface density profile 
S oc R~^, corresponding to a Mestel disk. In this frame¬ 
work, the effective viscosity of the disk follows from the re¬ 
quirement to maintain the accretion rate and the transport 
of angular momentum, t herefore prov i ding the amount of 
viscous heating (see e.g. Lodatdr2007t llnavoshi fc HaimanI 


l20l4[Latif fc SchleicheI^^2015a^ . We further distinguish be¬ 
tween scenarios where the rotation profile is determined by 
the self-gravity of the disk, as well as scenarios where grav¬ 
ity is dominated by the central source. 

If rotation is dominated by the central source, it fol¬ 
lows that the disk becomes gravitationally stable on smaller 
scales due to the relation S oc R~^. In particular in the 
regime of high accretion rates, one may however expect 
that self-gravity contrib utes signficantlv to the transp ort of 
angular momentum fe.g. lBegelman fc Shlosmai]l2009ll . We 
therefore consider also a self-regulated disk model, in which 
the disk viscosity is due to gravitational instabilities, there¬ 
fore requiring a scenario with the Toomre 1. In these 
cases, we obtain a steeper profile for the surface density 
of the disk in a marginally stable state. We refer here to 
iLodatol (1200711 for a review on the properties of such disks. 

The resulting evolution has been explored for a 


range of stellar masses and 
ering both conventional Poo. 

accretion rates, consid- 
III star formation fe.g. 

lAbel et al.ll2002tlBromm & Loebll2003tlYoshida et al.ll2008r) 

as well as the formation of supermassive primordial stars 

fe.g. Lodato & Nataraian 2006 

: Resmn & Haehneltl 120091: 

iLatif et al. I2013bl: iFerrara et al. 

[2014). As chemical initial 


conditions, we have both explored an initiall y molecular 
gas, a s expected for instanc e in minihalos (e.g. lAbel et al.l 
I 2 OO 2 I : iBromm &: LoebI 120031 : lYoshida et al.l l2008ll . as well 
as an initially atomic gas. In more massive haloes, the 
chemical evolution depends on the ambient radiation back¬ 
ground, and is typically molecular for moderate values and 
atomic for very strong values of the radiation background 
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(e.g. Omukail 200 lt Greif et alJ 200 ^ IShang et al.l[2ninbl : 


IScMeicher et alJl201C)HLatif et al.ll201lL l2015a l. 

In case of the initially molecular gas, we find that vis¬ 
cous heating plays a crucial role in the presence of a central 
source and accretion rates of ~ 10“^ Mq yr“^, both for 
the generic and the self-regulated disk model. For a cen¬ 
tral star with about 10 Mq, the viscous heating leads to 
the collisional dissociation of molecular hydrogen on scales 
of ^ 10 AU, and the effect is further enhanced for super- 
massive stars with ~ 10^ Mq, where the molecular gas 
is dissociated already on scales of ^ 1000 AU. The latter 
thus strongly contributes to the stabilization of the cen¬ 
tral region. Even if fragmentation occurs, the resulting 
clumps will be more massive, and are expected to rapidly 
migrate towards the central r egion ijlnavoshi fc HaimanI 
I20l4 iLatif fc Schleich^ l2015bll . The transition towards 
an atomic gas also increases the opacity, which is typically 
of order 1 cm^ g“^ at a temperature of 10^ K. The latter 
can have a relevant impact on the structure of the disk, and 
increase the overall stability. However, fragmentation may 
be more efficient on larger scales where the gas is molec¬ 
ular, and also the migration time is increasing with the 
distance from the central object. It is thus certainly vi¬ 
able that fragmentation occurs on some scale, which could 
potentially contribute to t he formation of a starburst ring 
around the central object (|Lodato fc Nataraianll20d^ . In 
such a case, the further accretion may be regulated by the 
interplay of gravitational i nstabilities with stell ar feedback 
(|Kawakatu fc Wadall2009l : IWutschik et al.ir2013ll . 

We further explored the case where the gas is initially 
atomic. In case of the generic disk model, it then remains 
atomic, as collisional dissociation of molecular hydrogen is 
efficient at the temperatu res of atomic hydrogen cooling 
(|Inavoshi fc Omukaill2012}l . Due to the high temperatures, 
the overall cooling rates are considerably enhanced, and 
the effect of viscous heating is therefore les s relev ant in this 
regime, as also described bv iFerrara et al.l l|2013f) for a hot 
metal-enriched gas. In the case of a self-gravitating disk 
model, we find that some H 2 may build up again towards 
the central region, due to the higher densities obtained in 
this regime. However, the characteristic temperatures are 
still above ~ 3000 K throughout the evolution. As the 
characteristic temperature for H 2 collision dissociation de¬ 
creases with increasing density, we find that the H 2 dis¬ 
sociates again after an initial build-up. Overall, the gas 
essentially remains atomic for most of the evolution. 

The potential impact of viscous heating is well-known 
for instance in the case of quasar accretion disks, where the 
resulting heat input leads to the presence of a hot gas over a 
large range of scales le.g. lGoodmanl[20n^lOoodman fc TanI 
l2004f) . In the context of primordial star formation, such 
effects have not been strongly considered. The one-zone 
models exploring the chemical evolution have typically as¬ 
sumed a free-fall collapse (e.g. lOmukai [200 It lOmukai et al.l 
[ 200 ^ [?ii over fc Ahell[200a [CH oveill2015h . while simulations 
were usually r estricted to t he ear ly stages of disk formation. 
For instance, IClark et al.l (1201111 were following the evolu- 
tion until ^ 1 _Mq of mass had turned into stars, while 
iGreif et'ffil l|201lh stopped the evolution when a 10 Mq 
star was formed, still with a moderate accretion rate of 
10“^ Mq. In a follow-up investigation at higher resolu¬ 
tion, they evolved the simulation until stellar masses of 

1 Mq, and noted that heating both due to the accre¬ 
tion onto the protostar, as well as due to shocks close to 


the spiral arm has occured (|Greif et al.l 1201^ . It is con¬ 
ceivable that these are the first signatures of viscous heat¬ 
ing, even though still at an early stage. In the context of 
supermassive stars, we note that the majority of runs al¬ 
ready started from an atomic gas, which typically remained 
atom i c during the furthe r e volution (e.g. Regan^_JIaehneIt 
20013: iLatif et al.l 2ni3aS iPrieto et al.l 20131 Regan et al. 
2014 : lBecerra et "l2015f) ~ The longer-term evolution in the 


presence of molecular coolin g ha s been explored in simula¬ 
tions by iLatif et ^ l|2014bfl and iLatif fc Volonter^ (l2015ll , 
showing that high accretion rates can be maintained, but 
without resolving the scales considered here. 


The results obtained here are particularly important for 
the formation of massive central object, as the enhanced 
temperature in the atomic gas increases the stability in the 
interior of the disk and helps to maintain a high accretion 
rate. The latter also alleviates the need of a strong external 
radiation field, typically expressed through a critical value 
Jcrit) as a means of providing an atomic gas. We expect 
these results to be particularly relevant in the case of mas¬ 
sive central objects of at least 10 Mq and high accretion 
rates of ^ 10“^ Mq yr“^. The nature of the phenomenon 
discussed here, and in particular the characteristic scaling 
of the viscous heating rate as R~^ leads us to the expec¬ 
tation that such a transition should generally occur, even 
though it may be shifted for instance if the surface density 
is enhanced and thus also the ability of the gas to cool. 
Even in the case of fragmentation and a more turbulent 
flow, we expect that gas velocities of the order of the Kep- 
lerian velocity will occur, leading to the formation of strong 
shocks and a transition towards the atomic cooling regime. 


Of course, the models here are still based on simplifying 
assumptions, such as the stationarity and axisymmetry of 
the disk, which can only hold at an approximate level. It 
is therefore important to investigate the evolution in such 
disks in 3D simulations following the interplay of chemistry, 
heating and cooling along with the gravitational dynam¬ 
ics to assess the impact on fragmentation. For such time- 
dependent models, the viscosity of the disk can no longer 
be obtained from the assumption of stationarity, but one 
may have to adopt explicit parametrizations of the viscosity 
i n self -gravitating disks as derived for instance by iRafikovI 
(|2015ll . Due to the timescales involved in the evolution of 
such disks, we expect that the latter can not be pursued in 
cosmological simulations, but that independent studies are 
needed to investigate the stability of self-gravitating pri¬ 
mordial disks with already a massive protostar. Depend¬ 
ing on the mass and accretion rate of the protostar, also 
radiation feedback may become relevant, even though it 
may predominantly affect the regions outs ide the disk (e.g. 
iHosokawa et al.l|201ll : iHirano et ffi] 12014!) . In the case of 
rapid accretion, it was found that the UV feedback from 
protostars is likely negligible, as the y are expected to have 
cool atmospheres like a red giant (IHosokawa et al.l 120131 : 
ISchleicher et al.l[2013ll . Previous investigations further have 
shown that a disk models with an optically thin cooling are 
therm ally unstable (jGannizzo fc Wheeler! 1984 iKato et Idl 
Il998ll . We note that the models pursued here are more 
complex, and in particular do not employ a constant a, as 
the disk viscosity is derived from Eq. [HI It is nevertheless 
conceivable that thermal instabilities may play a relevant 
role in their evolution, which needs to be explored in further 
detail. 
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Beyond the primordial case, the influence of dust 
cooling is potentially significant at high densities (e.g. 


Schneider et al.l 120031; lOmukai et al.l 20051: Schneider et al.l 

200(4 Omukai et al.l 

2008; Cazaux & Snaansl 20091: 

DoDcke et al.ll201lL 12012 

; Schneider et al.ll2012f). However. 


most of these investigations have employed one-zone mod¬ 
els assuming a free-fall collapse, while we have shown here 
that the heating mechanism is considerably different in a 
rotating disk around a central protostar. In this regime, 
the role of dust cooling requires further investigations, 
to understand the fragmentation behavior at the later 
stages of the evolution. The latter can be important 
both for the format i on of the first low -mass stars (e.g. 
ISchneider et al.lF2012t iKlessen et al.ll2012f) . but also for the 
formation scenarios of m assive black holes (|Omukai et al.l 
I 2 OO 8 I : l^rrara et al.ll201,‘l^ . Even if fragmentation occurs, a 
black hole may possibly form, as long as the cluster is fed 
via substantial inflows and gravitational instabili ties occur 
(jAlexander fc Nataraianl[20l4 iLatif et ^120159) . 
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Log Radius [AU] 




Fig. 1. Results for the generic disk models without central 
source and an initially molecular gas (Table [T]). Top panel: Gas 
temperature vs radius in case of a generic disk model without 
central source, assuming different accretion rates. Mid panel: 
Abundances of H, H 2 , H"*', H~ vs radius for the reference 
case NS3 (disk with no source, accretion rate 10“® Mq yr“^). 
Bottom panel: Heating and cooling contributions vs radius for 
the reference case NS3 (disk with no source, accretion rate 
10“^ Mq yr“^). We refer to Table [3] for the heating/cooling 
contributions defined in the legend. 
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Fig. 2. Results for the generic disk models with central 
source and an initially molecular gas (Table [21) • The species 
abundances and heating/cooling rates are given for the S2A3 
(generic disk model with central source of 100 Mq and accretion 
rate of 10“® Mq yr“^). We refer to the caption of Fig. [T] for the 
description of the panels. 
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Fig. 3. Results for the self-regulated disk models with an 
initially molecular gas (Table |4]|, plus species abundances and 
heating/cooling rates for the reference case SI (self-regulated 
disk with central source of 1 M© and an accretion rate of 
10“^ Mq yr“^). The top panel shows the gas surface density as 
a function of radius. We refer to the caption of Fig. [2] for the 
description of the other panels. 
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Fig. 4. Results for the generic disk model with an ini¬ 
tially atomic gas (Table [3. The species abundances and heat¬ 
ing/cooling rates are given for the HI (generic disk model with 
central source of 10 M© and an accretion rate of 10“^ Mq yr“^). 
We refer to the caption of Fig.[T]for the description of the panels. 
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Fig. 5. Results for the self-regulated disk models with an ini¬ 
tially atomic gas (Table [S|. The species abundances and heat¬ 
ing/cooling rates are given for the HSO (self-regulated disk with 
central source of 1 Mq and an accretion rate of 10~^ Mq yr“^). 
We refer to the caption of Fig.[3]for the description of the pan els. 
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